Fractional Quantum Hall Effect of Lattice Bosons Near Commensurate Flux 
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We study interacting bosons on a lattice in a magnetic field. When the number of flux quanta per plaquette 
is close to a rational fraction, the low-energy physics is mapped to a multispecies continuum model: bosons 
in the lowest Landau level where each boson is given an internal degree of freedom, or pseudospin. We find 
that the interaction potential between the bosons involves terms that do not conserve pseudospin, corresponding 
to umklapp processes, which in some cases can also be seen as BCS-type pairing terms. We argue that in 
experimentally realistic regimes for bosonic atoms in optical lattices with synthetic magnetic fields, these terms 
are crucial for determining the nature of allowed ground states. In particular, we show numerically that certain 
paired wave functions related to the Moore-Read Pfaffian state are stabilized by these terms, whereas certain 
other wave functions can be destabilized when umklapp processes become strong. 



Recent advances in the field of topological phases and their 
potential application in implementing an intrinsically fault- 
tolerant quantum computer ll] [2| have revitalized interest 
in fractional quantum Hall (FQH) states as the most promi- 
nent examples of topologically ordered phases of matter |3|. 
Even though it has only been observed in (fermionic) solid- 
statesystems, the FQH effect can also exist for bosons [4- 
l?). Promising candidates are systems of interacting ultra- 
cold atoms where the necessary magnetic fields are simu- 
lated by rapid rotation |8j |9J or by laser-induced synthetic 
gauge fields ifTOl . At low temperatures when the filling frac- 
tion ly (the ratio of the particle density n to the magnetic flux 
density n^) is sufficiently small, one can expect to observe 
bosonic versions of the FQH effect fSj. For example, the ex- 
act ground state of bosons with contact interaction at filling 
fraction i/ = 1/2 is the Laughlin state |i4jl11], while at = 1 
the ground state is in the same topological phase as the non- 
Abelian Moore-Read Pfaffian state lISHTIll. 

A major advantage of optical and atomic systems over con- 
ventional solid-state systems is the possibility of creating and 
controlling quasiparticle excitations more naturally and with 
higher precision (e.g., by shining focused laser beams on the 
atomic gas) ifTTI . A number of proposals suggest that the FQH 
regime for cold atoms can be most easily achieved using op- 
tical lattices flSWTTl. The question naturally occurs whether 
there is new physics that may arise for a system of interacting 
bosons in the FQH regime due to the effects of an underlying 
lattice. It has been shown that in the limit when the flux den- 
sity 710, or equivalently the number of flux quanta per lattice 
plaquette, is small, one can ignore the existence of the lattice 
and treat the system in the continuum limit | fSlH^. When 
is large, however, the presence of the lattice can potentially 
lead to new correlated states of matter that are absent in the 
continuum ||T9l42n . This is the limit we will focus on. 

The starting point for our analysis of the many-body 
physics in this problem is the observation that when is 
close to a rational fraction, the lowest energy bands in the Hof- 
stadter butterfly, a fractal structure realizing the single-particle 
energy spectrum of particles hopping on a lattice in a magnetic 
field |22J, are reminiscent of Landau levels in the continuum. 



This resemblance can be formalized by mapping the single- 
paiticle states of the system to a continuum model when the 
flux density is near simple rational fractions 1 19 1. 

The main result of this Letter is the following. In agree- 
ment with ||T9l , we find that for flux per plaquette close to 
a rational fraction, — p/q + e with p,q small integers, 
and e sufficiently small, one can map the system to an effec- 
tive continuum model with Landau levels and an added degree 
of freedom for the particles, a sub-band index or pseudospin, 
which can take q possible values. However, in addition to 
the density-density interactions between bosons of different 
pseudospin found in [19], we find anomalous "pairing" inter- 
actions that do not conserve the number of particles of each 
pseudospin species. We find that these pairing terms, which 
become increasingly strong as e is increased, are crucial in de- 
termining the possible ground states of the system for realistic 
values of the parameters of the problem. 

As a detailed example, we consider the most (experimen- 
tally) realistic case = 1/2 + e and study several effec- 
tive filling fractions i> = n/e. We find a new FQH state at 
D — 1, which does not exist without the pairing interactions 
but becomes stabilized by the increase in e and the concomi- 
tant increase in these interactions. This new state is related to 
two copies of the non-Abelian Moore-Read Pfaffian state (TT]; 
however, it is a topologically distinct phase of matter In 
contrast, we find that the pairing terms destabilize the states 
predicted at fillings D = 2/3 previously discussed by |T9l, 
D — 4/3 [20], and iy — 2 1,21 J . We present detailed numerical 
evidence for our conclusions and argue that experiments are 
most likely to be in a regime where these pairing terms are 
important. 

We consider bosons with onsite repulsive interaction, hop- 
ping on a two-dimensional square lattice, subject to a uni- 
form perpendicular effective magnetic field. This system is 
described by a modified Bose-Hubbard Hamiltonian 1231 . 

^ ^ --^ E {e''''4cj+h.c.) +Uj244c.c,, (1) 

<ij> i 

where cj and Ci are boson creation and annihilation operators 
on lattice site i, J is the hopping energy, and U is the strength 
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of the onsite interaction. Here 9ij = A ■ Al is the phase 
acquired by a particle hopping from site i to the neighboring 
site i with A being the vector potential, and we work in units 
where fi = 1, and the effective electric charge coupled to A 
is also set to unity. The kinetic term in the Hamiltonian then 
indicates that a particle hopping around a lattice plaquette ac- 
quires a phase of 2TTn^. 

We start by considering the kinetic term of the Hamiltonian 
only. This is the well-known single-particle Hofstadter prob- 
lem, which we review only briefly. We assume the lattice is in 
the X, y-plane, with lattice spacing set to unity for simplicity, 
and choose the Landau gauge so that A = (0, 27rn0 x, 0). The 
wave function becomes ip{x, y) — (j>{x)e^^'^ , where 4){x) sat- 
isfies Harper's equation (f){x + 1) + 4>{x — 1) + 2 cos(27rn0a; — 
k)(j){x) — {E/J)(f){x) and k is the momentum in the y- 
direction. Note that x and y are both integers. 

Consider the case of = e <C 1, where a continuum 
approximation of the discrete Harper's equation can be used 
for the low-energy eigenstates. In this limit, it is conve- 
nient to use a Wannier basis localized near minima of the 
cosine potential. These Wannier functions can be approxi- 
mated by harmonic oscillator (Landau level) solutions with 
oscillator length (magnetic length) = l/\/27re centered 
at Xk = k/{2TTe), i.e., (j)k{x) ^ exp(— 7re(a; — Xk)"^). The 
bandwidth of the lowest band arises from tunneling between 
adjacent minima of the potential and for small e it scales as 
^ e^'^l'^ where the constant C w 1.166 can be obtained by 
the WKB approximation |24|. Note that for small e, the band- 
width is much smaller than the band gap A — AJire, making 
this limit of the Hofstadter problem an example of an (almost) 
flat Chern band [25]. 

Now let us consider flux densities close to a rational frac- 
tion. For simplicity, we focus on = 1/2 + e, for 
which Harper's equation becomes (j){x + 1) + <j){x — 1) + 
2(— 1)"^ cos(27rea; — fc)(/)(a;) — {E/J)(t){x). This form sug- 
gests a Wannier solution analogous to the above case, but with 
a two-site form factor to account for the rapidly oscillating 
factor (—1)^. We propose the ansatz solution, 

~ (1 + A(-l)-+^)e-'^^(^-^'==)', (2) 

where we have defined Xks — {k — S7r)/(27re) and s = 0, 1 is 
a sub-band index. Thus for each momentum fc, there are two 
possible wave functions, which are spatially separated (due to 
the shift in the center of the oscillator) and also have their main 
weight on either the even or odd sites of the lattice. Choosing 
the value A = -\/2— l+7re(\/2— 2)/2 solves Harper's equation 
to order 0(e) and higher order terms can be added to A to 
satisfy the equation to still higher order. 

If we interpret the sub-band index s as a quantum number 
representing a new degree of freedom, the low-energy bands 
in the lattice at = 1/2 + e are equivalent to the energy 
bands of a two-species system at = e, which is the con- 
tinuum (Landau level) limit. Thus at rt^ = 1/2 + e we can 
understand e as the effective flux density giving rise to an ef- 
fective filling fraction defined as i> = n/e. Similarly, for gen- 
eral — p/q + e (with p and q coprime), q solutions can be 



found and the system can be treated as a q-species model with 
effective flux density e lfT9ll20l . With increasing q, the band- 
width increases as ^ e~^l'^'^ whereas the band gap, while 
remaining proportional to e, decreases with increasing q. 

In order to be in the FQH regime, it is necessary that the 
interaction energy be larger than the bandwidth (so that the in- 
teraction dominates over the kinetic energy). In addition, we 
would like the interaction to be smaller than the band gap so 
that all of the physics occurs within the lowest Landau band; 
however this requirement may not be crucial Il26l . Finally, 
the temperature must be less than the energy gap of the FQH 
state, which is typically set by the interaction energy (although 
it could also be set by the band gap if that is smaller). States 
competing with quantum Hall liquids include Bose-Einstein 
condensates: these describe the physics at — 1/2, for ex- 
ample 1271 128|. 

Given these restrictions, and given that experimentally ob- 
taining low temperatures will always be a challenge, it is clear 
that the FQH effect is most likely to be observed in the regime 
of intermediate e where, the band gap is not too small and the 
bandwidth is not too large. Indeed, it is perhaps optimal to 
work in a regime where bandwidth and band gap are compa- 
rable. One can simply look at the Hofstadter spectrum to see 
where these inequalities are best satisfied ll22ll . The most ex- 
perimentally favorable case occurs for = e <C 1. Here, e 
might be as large as 0.3 before the bandwidth is on the order of 
the band gap, and the band gap may be as large as J. This par- 
ticular case has been studied extensively previously ifTSlfTSl . 

The case of — 1/2 + e, which we focus on here, is also 
fairly favorable for the observation of FQH effect. The param- 
eter e can be as large as 0.1 before the band gap is on the order 
of the bandwidth, and the band gap may be as large as about 
0.3 J. While = 1/3 + e is still experimentally plausible, 
the cases of = p/q + e with q > 3 have extremely tiny 
band gaps and, hence, seem less accessible. We note that de- 
spite the fact that these inequalities of energy scales are harder 
to satisfy for — 1/2 + e than for = e, the former has 
richer physics associated with the new quantum number, the 
sub-band index introduced above. 

We now turn to consider the effect of the interaction term in 
the Hamiltonian Eq. ([TJ. Using any basis of (single-particle) 
states ipa{x,y) with corresponding creation and annihilation 
operators ip'l and i/'a, the interaction may be decomposed as, 

U=Y1 U-b-d ^l^UJd^ (3) 

ahcd 

where 

Uabcd = U ^ip*{x,y)ipl;{x,y)ipc{x,y)ipd{x,y). (4) 
x,y 

For n,p = e-^ 1, we use the basis for the lowest 
band, i.e., ijjkix^y) — (j)k{x)e^^y with (pk{x) the Gaus- 
sian form as described above (properly normalized). In this 
limit we may convert the sums into integrals, then, pro- 
jected to the lowest energy band, we obtain Ukik2k3ki = 
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FIG. 1: (a) The FQH gap E at effective filling fraction ^ = 1 as 
a function of e, where the flux density is — 1/2 + e. Data are 
shown for A'' = 8, 10 bosons. Increasing e increases the strength of 
the pairing terms of the Hamiltonian and stabilizes this state, (b) The 
FQH gap for i> = 2/3 (N = 10), 4/3 (iV = 12), and 2 (iV = 14) as 
a function of e. Increasing e decreases the FQH gap and destabilizes 
the corresponding states, (c) Overlap between the exact ground state 
of the system at effective filling fraction i> — 1 and the trial wave 
function Eq. |6]l vs. e. The overlap for e > 0.1 exceeds 95 %. 



ViC/e-^-<^=i('=*-^-^)'/(i6"^)4i+fc.-fc3-fc4. where the func- 
tion Sp is defined to be unity if the argument p is an integer 
multiple of 2tt and is zero otherwise. The Gaussian factor en- 
forces ki w kj so that total momentum fci + ^2 — ^3 — ^4 must 
be zero, not just mod 27r. This derived form of the interac- 
tion is precisely what we expect for continuum bosons in the 
lowest Landau band with short-range contact interaction f29\. 

For = 1/2 + e, we can use the basis ijjks{x,y) = 
<j>ks{x)e^''y with (pks{x) given by Eq. Projecting to the 
lowest Landau level, we correspondingly find, 

C^/ciSife2S2fe3'S3'=4S4 ~ 

_ 4 (fci-fcj^7r(3.i-Sj))^ 

where the matrix G results from summing over the discrete 
form factors in the expressions for the corresponding wave 
functions. Note that the functional form of the interaction 
Eq. (|5]l is identical to contact interactions for a continuum 
Landau level up to the band index dependent matrix G out 
front once we redefine the momentum as fc = fc — its. In 
terms of these new variables, the Gaussian enforces ki « kj 
which now allows ki + k2 — kj, — k^ — ±27r (allowed by 5) 
if Si + S2 — S3 — S4 = ±2. 

Given the precise resemblance to a two-species continuum 
model, we may employ exact diagonalization, which is well 
established as a numerical technique for the study of interact- 
ing particles in continuum Landau levels 1291 . We chose to 
formulate the effective problem obtained from our preceding 
analysis in a finite spherical geometry |30, 31 1, which elimi- 
nates edge effects and thus allows direct access to the physics 
of the bulk. This problem is defined entirely within the lowest 
energy band by the interaction, Eq. (|5]l, as the kinetic energy is 
assumed to have small bandwidth compared to the interaction 
scale U . 

Following |fT9l we switch to a new basis ipks — {4'ko + 



s i'4'ki)/V^, where s = ±1 (or 'up' (f) and 'down' (\.) in the 
current text). We refer to this new form of the sub-band index 
s as pseudospin. Using this basis, the nonzero elements of the 
transformed matrix G become G^^^^ = Gim = G^i^i = 
^Itlt ~ ^Ittl = ^-^tllt ~ 1- addition, we also find 
two extra "pairing" terms: G^^^ — G^^n — ne. These 
terms correspond to fci + A:2 = ^3 + ^4 + 27r, which resem- 
ble umklapp scattering processes, and do not conserve pseu- 
dospin (si + S2 ^ S3 + S4). These terms indicate that a 
pair of pseudospin ups (downs) can be annihilated while a pair 
of pseudospin downs (ups) are created thus suggesting that a 
BCS-type pairing |[32l[33l could occur between particles with 
the same pseudospin. While these pairing terms vanish in the 
limit of e — J- 0, as mentioned above, the experimentally rele- 
vant regime is likely to be at finite e where these terms will be 
important. 

The novel twist for lattice bosons near = 1/2 is the 
emergence of an umklapp scattering term between the two 
emergent species. Consequently, the total pseudospin is not 
conserved, and we need to take into account the full Hilbert 
space containing all possible distributions of particles into 
the two sub-bands. We examine the spectrum of the Hamil- 
tonian Eq. (j5]l for the occurrence of incompressible ground 
states that are characterized by translational invariance (angu- 
lar momentum L = on the sphere) and a finite FQH gap. 
Our search yields four candidates at effective filling fractions 

= n/e = 2/3, 1,4/3 and 2 where = l/2 + e[,34J. 

For the effective filling fraction i> = 1, we find that the en- 
ergy gap between the ground state and the first excited state 
rapidly opens up as one increases e [see Fig.[T|a)]. This indi- 
cates that the umklapp pairing terms, G^^^^ and G^^^^, which 
are the only terms in Eq. (j5]l that change with e, are respon- 
sible for producing an incompressible state at this filling fac- 
tor. On the contrary, the energy gaps at effective filling frac- 
tions i> = 4/3 and v — 2 close as e increases, indicating that 
the pairing terms destabilize the corresponding incompress- 
ible states [see Fig.[TJb)]. Atv — 2/3, as pointed out in 1 19|, 
the 221 state is an exact ground state, and this remains true 
even in the presence of the pairing terms. However, as e in- 
creases, this gap also closes, as excited states are sensitive to 
the pairing terms [see Fig.[TJb)]. 

Interestingly, the energy gaps of the i> = 2/3 and v = 2 
states have a very similar magnitude and dependency on e, 
which may be unexpected in the current formalism. This, 
however, is a natural conclusion in the composite fermion 
(CF) theory for the lattice, which explains both these states 
by the same energy gap in the CF spectrum ||2T]. Further- 
more, our expansion of the effective model to linear order in 
e predicts that the gap of these states closes near the value, 
6 = 1/6, predicted by CF theoiy ED- The case of V = 4/3 
is clearly very different. Here we find that the ground state at 
small e has significant overlap of > 0.7 with the non-Abelian 
spin singlet (NASS) state 1351 for = 12 particles. However, 
this state is very fragile and we cannot ascertain that it persists 
in the thermodynamic limit. 
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As was mentioned above, the presence of the umklapp pair- 
ing terms in the Hamihonian, G^^^ and G^^n, suggest a 
BCS-type pairing between particles of the same pseudospin. 
To further investigate the effect of these terms on the nature 
of the ground state at i> = 1 we propose an explicit trial wave 
function. We use the common conventions for studying the 
FQH effect: we adopt the symmetric gauge and write coordi- 
nates of particles on the plane in dimensionless complex form 
z ^ {x + iy)/lct- Our trial wave function is given in the co- 
ordinates of bosons of 'up' type and bosons of 'down' 
type (with both A'^ and A'^ assumed to be even): 

Nt- 

i<j=i i<j=i j=i 

As written, the up and down bosons are assumed distinguish- 
able. Expanding this wave function in the original bosonic 
basis, one obtains an expression that is fully symmetric in all 
coordinates. 

Here, the two Pfaffian factors (Pf) are antisymmetrized 
sums over pairs of particles with the same pseudospin 
Pf(j-^) = -^(i73ijr^i7 • • ■ ) with A the antisymmetriz- 
ing operator. This Pfaffian form is precisely the real space 
form of a BCS pairing wave function, which indicates that 
particles with the same pseudospin form pairs. Without the 
Jastrow factors this type of pairing is analogous to a ^He-A 
phase with a d-vector in the x, y-plane |32|. As in the case 
of other paired Hall states, the topological properties are only 
trivially altered by restoring the Jastrow factors ||32l . 

Note that the wave function, Eq. (|6]l, is the exact ground 
state of H2-3, the sum of a three-body delta-function inter- 
action for particles with the same pseudospin and a two-body 
delta-function interaction between particles of opposite pseu- 
dospin. 

To check the validity of this trial wave function (obtained 
as the ground state of H2-3 on the sphere), we calculate its 
overlap with projection of the ground state of Eq. (|5]l onto the 
sector with — N^. As e is increased and the gap opens 
up, we find increasing overlap between our trial state and the 
exact ground state [see Fig. [TJc)]. The overlap for e = 0.1 
is above 95 % for iV = 12 particles, which is an excellent 
indicator of the validity of our proposed wave function. Note 
that, although outside the regime of validity for our model, at 
e ~ 0.16, Eq. (j6]l is nearly an exact ground state of the two- 
body interaction, Eq. (jsj), to an accuracy of about 10^^. We 
also find that the inverse 1 /z dependence of the paired wave 
function in Eq. ^ is optimal, as introducing variational pa- 
rameters to change its shape [33] does not increase the overlap 
significantly. 

We have also studied the quasihole spectrum of Eq. (j6]l in 
the presence of additional flux. For the model H2-3 Hamil- 



tonian, the quasihole spectrum is precisely that of two decou- 
pled Moore-Read layers — the quasiholes of each layer cor- 
responding to the so-called half-quantum vortices of '^He-A. 
However, for our Hamihonian of interest, Eq. (j5]), the umk- 
lapp pairing terms lock the direction of the d-vector, thus re- 
quiring that the quasiholes pair between layers, confining the 
half quantum vortices and leaving the system with effectively 
Abelian excitations. To establish with clarity that this is the 
correct physics we have been able to predict the entire low- 
energy quasihole spectrum of the ■^He-A model, using a gen- 
eralization of the approach introduced in [36 J , which precisely 
matches the low-lying spectrum of the microscopic Hamil- 
tonian, Eq. (|5]l, for every case we could numerically access. 
These results will be presented elsewhere. 

Signatures for our proposed state can be derived from a 
range of experimental probes for the detection of quantum 
Hall states in cold gases, such as measurements of ground- 
state incompressibility [3?J, noise correlations [20, 38J, and 
possibly a direct measurement of quasihole statistics ifTTIl . 

The methods described here can be generalized to flux den- 
sity = p/q + e although, as discussed above, FQH effect 
with larger q is likely to be harder to realize in experiments. In 
this case there would be a sub-band index s — 0,1 . . . {q — 1) 
and the umklapp terms of the interaction would allow noncon- 
servation of this sub-band index via si + S2 — S3 ^ S4 = 
mod q, which could lead to new pairing terms and possibly 
new physics. 

To summarize, we have shown that anomalous pairing 
(umklapp) interaction terms are crucial to the physics of FQH 
effect for interacting bosons on a lattice at flux density — 
1/2 + e. We find that the pairing terms greatly modify the 
ground state at various effective filling fractions i> = n/e. 
At i) — 1, we demonstrate that these terms stabilize a new 
paired FQH state, which is effectively two coupled copies of 
the Moore-Read Pfaffian state. At i> = 2/3, 4/3 and 2, we find 
that the incompressible states are destabilized by the pairing 
terms. 
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